Metastable superfluidity of repulsive fermionic atoms in optical lattices 
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In the fermionic Hubbard model, doubly occupied states have an exponentially large lifetime for 
strong repulsive interactions U . We show that this property can be used to prepare a metastable s- 
wave superfluid state for fermionic atoms in optical lattices described by a large-f/ Hubbard model. 
When an initial band-insulating state is expanded, the doubly occupied sites Bose condense. A 
mapping to the ferromagnetic Heisenberg model in an external field aUows for a reliable solution of 
the problem. Nearest-neighbor repulsion and pair hopping are important in stabilizing superfiuidity. 

PACS numbers: 03.75.Ss,71.10.Fd,67.25.D- 



Trapped cold atoms open the possibility to realize new 
quantum states of matter and to control them with an un- 
precedented precision. An especially exciting perspective 
is the possibility to study interacting quantum systems 
out of equilibrium. The high tunability in combination 
with the slow dynamics of cold atoms allows to investi- 
gate time-dependent processes, for example the quench 
from a superfluid to a Mott insulating state [1] . 

Thermal equilibrium is usually dominated by low- 
energy states of the system, while out of equilibrium also 
high-energy states can become important. In the con- 
tinuum, high-energy states typically decay rapidly in the 
presence of interactions: High energy implies that the 
available phase space for inelastic scattering is large. In 
contrast, for lattice systems where the kinetic energy of 
a single particle cannot exceed its bandwidth D, a state 
with high energy, E ^ D (e.g. a doubly occupied site in a 
strongly repulsive Hubbard model), cannot easily decay. 
This is a consequence of energy conservation: To dissi- 
pate the huge energy E, a complex many-particle scat- 
tering process is needed, with at least n > E/D partici- 
pating particles. For local two-particle interactions such 
processes are expected to be exponentially suppressed for 
large n (see below) . This effect has been directly observed 
in measurements of the lifetime of doubly occupied lat- 
tices sites for bosonic ®^Rb atoms in an optical lattice [2] : 
Starting from a dense cloud of atoms with many doubly 
occupied sites, the strength of the trapping potential was 
reduced in one direction, allowing the cloud to expand. 
Subsequently, many long-lived double occupancies were 
detected, with a lifetime exceeding their inverse tunnel- 
ing rate by more than two orders of magnitude. 

The large lifetime of doubly occupied lattice sites im- 
plies that one can easily create new metastable states 
of matter. Indeed, numerical simulations by KoUath et 
al. [3] show that metastable states form in the one- 
dimensional bosonic Hubbard model for strong repulsion. 

An obvious question is whether the doubly occupied 
sites will Bose condense. For a bosonic Hubbard model, 
this question was investigated by Petrosyan et al. [4] , but 
the authors found that instead the system will phase- 
separate: Due to nearest-neighbor attractive interac- 
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FIG. 1: Schematic plot of the momentum distribution of 
fermion pairs [7]. For attractive interactions, the Cooper 
pairs condense at momentum (and corresponding recipro- 
cal lattice vectors). In contrast, the metastable superconduc- 
tivity of the repulsive Hubbard model arises at momentum 
(±7r, ±7r, ±7r). 



tions, doubly occupied sites will stick together instead 
of forming a low-density superfluid. In this paper we 
will prove that for fermions, in contrast, a Bose conden- 
sate of spin singlets with s-wave symmetry will form. 
Interestingly, the many-particle wavefunction of the rel- 
evant homogeneous metastable superfluid state can be 
constructed in a controlled way. It has been known for 
a long time [5, 6] that a hidden SU(2) symmetry of the 
charge sector (called SUc (2) in the following) of the Hub- 
bard model can be employed to build wavefunctions with 
off-diagonal long-range order (states with so-called "t? 
pairing" [5]). We shall show that these states can eas- 
ily be realized just by expanding an atomic cloud in an 
optical lattice slowly compared to typical collision times 
but rapidly compared to the exponentially large lifetime 
of the doubly occupied states. 

The condensation of doubly occupied sites can be 
detected by measuring the momentum distribution of 
fermion pairs [7]. The repulsively bound doubly occu- 
pied sites of the repulsive Hubbard model hop from site 
to site via virtual low-energy states. Therefore the sign of 
their effective hopping amplitude is reversed compared to 
bound pairs in the attractive Hubbard model. This im- 
plies that the condensation occurs at momentum (tt, tt, tt) 
[5, 6] rather than zero, allowing for an unambiguous de- 
tection of this state, see Fig. 1. 
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Setup. We consider the fermionic Hubbard model 

{ij),(T i i 

on a cubic lattice with a harmonic trapping potential 
of strength Vt- Here J is the tunneling rate between 
neighboring sites of the optical lattice, U > 12J = D 
is a strong repulsive interaction, and riia- = c\^Cia, rii = 
nil + "^«T- '^^^ lattice distance o is set to unity. 

As argued above, the total number of doubly occu- 
pied sites, = ?iit?iij.' extremely long life- 
time, due to the difficulty in loosing the large energy U. 
This can formally be seen using a well-known unitary 
transformation [8, 9], H H = e^^Ue-'^ , Nj. ^ Nd = 
e^^ Nde~^^ , called Schrieffer- Wolff transformation. For a 
given arbitrary order n, one can explicitly construct [9] a 
unitary operator e*'^" , such that the commutator [Nd, Ti] 
vanishes exactly up to terms of order l/{7". This proves 
that, in the limit of large U, the lifetime of doubly 
occupied sites grows faster than any power of U. The 
underlying physical reason, the energy bottleneck, has 
been described in the introduction. We therefore expect 
that, for Vt = 0, Td is exponentially large in U/D. 

We consider an initial situation where the atoms are 
densely packed, with two atoms per site in the center 
of the trap (i.e. a band insulator state), and investigate 
the evolution of the system upon reducing the strength 
Vt (i.e. curvature) of the trapping potential [10]. The 
initial system is in thermal equilibrium, and we assume 
vanishing entropy for simplicity (all of the following ar- 
guments remain valid as long as the entropy per particle 
remains small compared to unity). The radius of the 

1/3 

atomic cloud is ~ A^^ . To avoid a decay of the dou- 
bly occupied states by a conversion of interaction energy 
into potential energy, the slope of the trapping potential 
at the edge of the cloud has to be small compared to 
U, 2Vtrd <C U. Taking into accoimt the Mott-insulating 
shell forming around the band-insulating core [11, 12, 14], 
one obtains from this condition the ratio of the numbers 

1 /3 

of singly and doubly occupied sites, Ni/Nd » l/^d • 
Nevertheless, the ratio N\/Nd can be made sufficiently 
small, such that singly occupied states can be neglected. 
Note that this is not required to obtain Bose condensa- 
tion of double occupancies, but simplifies the theoretical 
analysis considerably. 

Effective m,odel. Neglecting singly occupied sites, the 
effective Hamiltonian after the Schrieffer- Wolff transfor- 
mation [8, 9] reads [15] (up to constant contributions) 




Wt^r^n^tn^i. (2) 

i 

The first term describes the hopping of doubly occupied 
sites, the second an effective interaction. In the pres- 



ence of singly occupied sites, the leading correction to 
(2) arises [9] from Jj2{ij)a''^h-cr'^la'^jcrnj,-a, which de- 
scribes an exchange of a doubly and a singly occupied 
site. This term can be neglected when the local density 
of single occupancies is smaller than J/U. While this is 
not the case at the border of the atomic cloud in its ini- 
tial configuration, it turns out to be valid in the scaling 
limit discussed below, as single occupancies are efficiently 
diluted when the trapping potential gets weaker. 

It is useful to rewrite (2) in two different ways. First, 
one can identify the doubly occupied states with a boson 
dl = cj-|.ctj^ such that (up to a constant) 

J2 

{ij) ij i 

with Udi = d\di. Here Va = oo implements a hard-core 
constraint, and Vij = —J'^/U describes an attraction for 
nearest neighbors i and j. Second, one can map the hard- 
core bosons to spins [4]. Starting from (2), this can be 
done by performing a particle-hole transformation for the 
down-spins only, (^^^ cjp cl^ — > (— l)'Cj|. This maps 
an empty site to a spin down and a doubly occupied site 
to a spin up. A finite magnetization in the xy plane de- 
scribes Bose condensation of pairs of fermions, see below. 
Identifying Sj = 5 X^a/s ^L"'"/?^*/? softer this transforma- 
tion, one obtains a ferromagnetic Heisenberg model in a 
magnetic field: 

n = ~^Si-Sj + 2Vt^rtSt. (4) 

(ij) i 

The SU(2) symmetry of the first term in (4) is a direct 
consequence of the SUc(2) symmetry of the charge sector 
of the underlying Hubbard model [5, 6]: For Vt = and 
a chemical potential /x = U/2, 7i (1) commutes with all 
three components of the particle-hole transformed op- 
erators Si defined above - in the original variables, 
these are (r; + r;t)/2, {rj -r]'')/{2i), and EiK-l)/2 with 

For Vf = 0, the exact ground state of (4) for fixed 
particle density rid is a ferromagnetic state, 

I*) = e-*^^*'^?| TTT •••) = e-'^^'(-^^*('=n'^*i+''-^-^|0) (5) 

where cos 6 = 1 — 2nd fixes the magnetization Sz =nd — \ 
in z direction. The rotational symmetry around the z 
axis is spontaneously broken, which implies off-diagonal 
long-range order, (c||c||) = {d\) = sing, with mo- 
mentum (tt, 7r,7r). The superfluid fraction, defined as 
\{d\)\'^ /rid-, is given exactly by (1 - rid) [16]. 

The initial thermally equilibrated state, described 
above, is not superfluid. A state with a finite expectation 
value of (cj|c||) can, however, easily be generated by re- 
ducing the trapping potential Vt adiabatically, i.e., much 



FIG. 2: Density profiles of double occupancies in the trap for 
different values of the trapping potential Vt , expressed as at = 
VtN^^U/ , a) for vanishing and b) for finite next-neighbor 
repulsion and pair hopping, 7 = {6U' + 3Jp /2)U N'^^^ / J'^ , see 
text. Upon reducing Vt, the cloud expands and, simultane- 
ously, the double occupancies condense, |(c|c|)p/nd = l — rid- 
In b), the atomic cloud spreads out further and superfluidity 
is enhanced. 



slower than the typical time scales of order [// aris- 
ing from the dynamics of the effective Hamihonian (2). 
If the whole experiment is furthermore performed on a 
time scale short compared to the (exponentially large) 
life time of doubly occupied sites, the system remains 
in the ground state of the effective Hamiltonian (2) [or, 
equivalently, (4)] with the number of doubly occupied 
sites (the total magnetization) kept fixed. To calculate 
how a finite superfluid fraction arises when Vt is low- 
ered for fixed Nd, we employ a mean- field approximation. 
We approximate the ground-state wavefunction of (4) by 
^P) = Ililni), where \ni) describes a spin i polarized in 
the +'hi direction. Here fii — h{vi) is (in the limit of 
large N and small Vt) a unit vector smoothly varying as 
a function of r, which minimizes 



L.(yn{v)f + Vtv\n^v) + l) 



(6) 



with the constraint / d^r[n^(r) -f 1] = '^Nj,. The em- 
ployed continuum limit becomes formally exact for a 
large number of atoms (typical optical-lattice experi- 
ments use 10"* to 10^ atoms). In this "scaling" limit it 



is convenient to measure distances in units of N, 



1/3 



and 



all physical properties only depend on the dimensionless 
quantity at — VtN^^^U / J'^ . The initial state then corre- 
sponds to at = cx). 

Results. When the confining potential Vt gets weaker, 
the cloud of atom pairs expands (Fig. 2a), and simul- 
taneously a condensate of fcrmionic pairs first emerges 
at the boundary of the band-insulating state, i.e. at 
the domain wall separating the spin-up and spin-down 
phases (Fig. 3a). When at becomes of order 1 and 
smaller, the maximum of |(cjc|)| moves to the center 
of the trap, and for at — > the condensate fraction 
shown in Fig. 3b approaches unity. The minimization of 
(6) then becomes equivalent to solving the Schrodinger 




FIG. 3: a) Square of condensate amplitude, |(c|^cj)p, as func- 
tion of the radius and the strength of the trapping potential 
Oft, for next-neighbor repulsion and pair hopping 7 = 25. 
b) Superfluid fraction as function of the trapping poten- 
tial, plotted as a~[^^^ . c) Projected momentum distribution 
nd{k±) = / dkz{d[,dk) , for different at as function of the dis- 
tance to (tt, tt). The condensate leads to a characteristic peak, 
while the non-condensed part gives rise to a uniform back- 
ground (not shown). The oscillations for large a are signa- 
tures of the superfluid shell in the early stage of the expansion. 



equation of non-interacting bosons in a trap using fi « 
(M/(r),0,-l + vE'(r)V2). 

Is the system a true superfluid? A superfluid is very 
different from a ferromagnet, as the excitation spectrum 
of the former is linear in momentum while it is quadratic 
in the latter. As the bosonic Hamiltonian (2) is equiv- 
alent to a ferromagnetic model, we conclude that, for 
14 = 0, the doubly occupied states do not form a su- 
perfluid with a finite phase stiffness, but only a Bose- 
Einstein condensate as for non-interacting bosons: At 
low energy and density, the scattering length of the 
bosons vanishes, as the hard-core repulsion is exactly 
balanced by the short-range attraction. Furthermore, the 
SUc(2) symmetry of the Hubbard model implies that the 
energy difference per volume of a phase-separated state 
(where doubly occupied sites stick together, i.e. where S 
points only up or down) and a superfluid state vanishes 
in the thermodynamic limit of the uniform system. 

SUc(2) symmetry breaking will therefore either sta- 
bilize the superfluid state or lead to phase separation. 
Let us list possible corrections to the Hubbard model (1) 
which break the SUc(2) charge symmetry (at Vt = 0), 
but preserve the SU(2) spin symmetry. As the chemi- 
cal potential is fixed by the particle number, the most 
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important contributions are two-site terms. The possi- 
ble two-site terms are longer-range tunneling, assisted 
tunneling, pair tunneling, density-density and spin-spin 
interactions. For an optical lattice, where the lat- 
tice potential in x, y, and z directions is independent, 
^{^) = the leading longer-range tunneling 

term is to the second-neighbor site in x, y and z di- 
rection. Its strength can be estimated as J' ~ J^/A 
where A measures the gap to the next Bloch band of 
the lattice. Assisted next-neighbor tunneling of the form 
—Jcl^Cjc,{ni^-cr + nj-a — 1) arises from the interaction 
correction to the local Wannier wavefunction of a fermion 
and is hence given by J JU / A. Finally, both the next- 
neighbor density-density interaction, U'{ni — l){nj — 1), 
and a next-neighbor pair-hopping term JpC^i-^c^i]^CjiCj'\ are 
given by matrix elements of Wannier states on adjacent 
sites, U' « Jp/'^ ~ U{J/A)^. A next-neighbor spin-spin 
interaction is of similar size, but unimportant for the dy- 
namics of doubly occupied and empty sites. One may 
also consider three-site terms, but those are easily seen 
to be subleading. 

After the Schricffcr- Wolff transformation, the leading 
correction to the dynamics of double occupancies for 
J <C ?7 ^ A arises from the U' repulsion term and the 
pair hopping Jp, as the tunneling terms only enter in sec- 
ond order via an intermediate singly occupied state. In 
fact, the contributions of J cancel, and the effect of J' 
is J''^/U - J^/{A'^U), such that the effects of U' and 
Jp are larger by a factor of (U/J)'^ » 1. Hence, in the 
bosonic language we are left with 

An = J2 4t^'K,i - ^)Ki - ^) + Jpdldj (7) 

iv) 

and both terms stabilize superfluidity. Translating into 
spins and using the same variational Ansatz as above, 
one obtains the leading correction to (6) in d = 3: 
AE = (6C/' + 3Jp/2) / d^rfil where we used = 1. 
It is convenient to parameterize the strength of the 
additional interactions by the dimensionlcss parameter 
7 = (QU' + ■iJp/2)UNy^/j'^. For a typical experimental 

setup 7 - {U/AfNl'^ will be quite large, 10 < 7 < 
1000. As shown in Figs. 2b and 3b, a large 7 leads to 
the expected expansion of the cloud and therefore to an 
enhancement of the superfluid fraction. 

From our variational wavefunction, one can calculate 
the momentum distribution of the fermion pairs, {d\,dk)- 
While the non-condensed fraction gives only a smooth 
background signal, the condensate gives rise to sharp 
peaks centered at (±7r, ±7r, ±7r), see Fig. 1 and Fig. 3c. 

Conclusion. We have shown that, within the strongly 
repulsive Hubbard model, one can realize a metastable s- 
wave superfluid by expanding a band-insulating ground 
state. One experimental problem may be the preparation 



of the initial state, as in present optical-lattice experi- 
ments [13, 14] the entropy is not small, implying a size- 
able fraction of singly occupied sites [14] even for large 
Vi . Fortimately, the condition for the onset of Bose con- 
densation is not very strict. For non-interacting bosons, 
the entropy per particle has to be smaller than 3.6^3. 
The corresponding entropy per fermion of l.Sfcs can be 
reached by cooling non-interacting fermions down to 0.22 
Tp; considerably lower temperatures are nowadays ob- 
tained routinely [14]. It is, however, presently not clear 
how this entropy is distributed between singly and dou- 
bly occupied sites. Furthermore, it may not be simple 
to keep quasi-adiabatic conditions when the cloud is ex- 
panded [13, 14], as extrinsic heating processes limit the 
total time in which experiments are performed. While 
a quantitative estimate of these corrections is difficult, 
we are optimistic that, with present-day technology and 
suitably chosen experimental conditions, an exotic finite- 
momentum s-wave condensate can be realized and de- 
tected for strongly repulsive fermions in optical lattices. 
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